Wave chaos in rapidly rotating stars 
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Effects of rapid stellar rotation on acoustic oscillation modes are poorly understood. We study 
the dynamics of acoustic rays in rotating polytropic stars and show using quantum chaos concepts 
that the eigenfrequency spectrum is a superposition of regular frequency patterns and an irregular 
frequency subset respectively associated with near-integrable and chaotic phase space regions. This 
opens new perspectives for rapidly rotating star seismology and also provides a new and potentially 
observable manifestation of wave chaos in a large scale natural system. 
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Since helioseismology revolutionized our knowledge of 
the Sun's interior, many efforts including space missions 
(MOST, COROT and KEPLER) are undertaken to de- 
tect oscillation frequencies in a large variety of stars [HQ. 
But to access the information contained in these data, 
the observed frequencies must be first associated with 
the right stellar oscillation modes. This crucial identifi- 
cation process requires a full understanding of the prop- 
erties of the oscillation spectrum and, for slowly rotating 
stars like the Sun, the asymptotic theory of high fre- 
quency acoustic modes provided such an understanding 
Both the approximate treatment of the centrifugal 
distortion [3j and the lack of asymptotic theory have so 
far hindered reliable identifications in rapidly rotating 
pulsators. This long-standing problem mainly concerns 
massive and intermediate-mass stars [3] like the 8 Scuti 
star Altai'r whose surface oblateness has been measured 
by interferometry Accurate computations of acous- 
tic modes fully taking into account the effect of rotation 
on the oscillations have only recently been performed for 
polytropic models of rotating stars [6| . Here we construct 
the dynamics of acoustic rays to understand the proper- 
ties of the frequency spectrum. 

The acoustic ray model is analogous to the geometrical 
optics limit of electromagnetic waves or the classical limit 
of quantum mechanics. The construction of eigenmodes 
from stellar acoustic rays has been already considered in 
the integrable case of a non-rotating spherically symmet- 
ric star [3] • However when the ray dynamics is no longer 
integrable, the problem is known to become of deeply 
different nature. This issue has been mostly investigated 
by the quantum chaos community in the context of the 
classical limit of quantum systems Q and the developed 
concepts have been applied to other wave phenomena 
such as those observed in e.g. microwave resonators 
lasing cavities [lo| . quartz blocks [ll[ and underwater 
waves 12]. The potential interest for stellar seismology 
has been suggested [l3| but not yet demonstrated. 

Our star model is a self-gravitating uniformly rotating 
monoatomic perfect gas (r = 5/3) where pressure and 



density follow a polytropic relation P e oc pl +1 ^ N with 
N = 3. Neglecting the Coriolis force and the gravita- 
tional potential perturbations, small amplitude adiabatic 
perturbations around this equilibrium verify: 
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where the density p, pressure P and velocity v describe 
the perturbation while c s is the sound velocity and g e 
is the effective gravity resulting from the gravitational 
and centrifugal potentials. As other quantities charac- 
terizing the star model, c s and g e vary in the merid- 
ian plane of the rotating star. Neglecting gravity waves, 
these equations can be reduced to the form, (w 2 — lj 2 )$ — 
c 2 A$ = where $ = P/c? s is related to the time- 
harmonic pressure perturbation P = 5R{Pexp(— iujt)} 
and uj c = v /(15/64)(g e /c s ) 2 + (3/8)V-g e is the cut-off 
frequency whose sharp increase in the outermost layers 
of the star provokes the back reflection of acoustic waves. 
The WKB approximation then leads to the eikonal equa- 
2 . The acoustic ray is the trajec- 
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tory tangent to the wave vector k at the point x and 
its evolution can be described by Hamilton's equations, 
H = c 2 k 2 + lo"1 being the Hamiltonian pj. Rays head- 
ing towards the star center tends to be refracted by in- 
creasing sound velocity while close to the surface non- 
specular reflection takes place at lo c = lu. As rotation 
increases, isocontours of c s and uj c are distorted together 
with the star surface. In the following, we restrict our- 
selves to axisymmetric modes L z = 0, thus reducing the 
phase space to four dimensions. 

The acoustic ray dynamics has been investigated by in- 
tegrating numerically the Hamilton's equations and the 
resulting dynamics is visualized using the standard tool 
of Poincare surface of section (PSS). As rays do not reach 
the star boundary, the PSS is defined by r p (0) — r s (9) — d 
where the distance d — 0.08r s (7r/2) from the star surface 
r s (6) has been chosen such that all but a few whisper- 
ing gallery rays cross the PSS (only outgoing rays are 
taken). The two coordinates of the PSS are 9, the co- 
latitude, and kg/to, kg being the angular component of k 
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FIG. 1: (Color online) PSS and typical acoustic rays at a rotation equal to 59 percent of the Keplerian limit. A whispering 
gallery ray (green/light grey), an island ray (blue/dark grey) and a chaotic ray (red/grey) are shown on the physical space and 
on the PSS (diamonds in the center figure). The central orbit of the island is also shown (yellow between points a and 6). kg/u 
is in unit of ^/rf (0)/GM with M the stellar mass and G the gravitational constant. 



in the natural basis associated with the coordinate sys- 
tem (( = r s (9) — r, 6). We use the scaled variable kg/u) 
as, in the limit uj ^> oj c , the ray dynamics becomes inde- 
pendent of the frequency away from the reflection points. 
We found that increasing the stellar rotation leads to a 
soft transition from integrability to chaos analogous to 
the one described by the KAM theorem. As illustrated 
in Fig. 1 for a given rotation rate, the phase space shows 
a mixed structure where chaotic regions coexist with a 
whispering gallery region close to the boundary and reg- 
ular islands around stable periodic orbits. As rotation 
increases, both the chaotic region and the central island 
chain get larger. A crucial feature of the dynamics is that 
each region is dynamically isolated from the other by in- 
variant tori which prevent communication between them. 
Such a situation has been found several times in the do- 
main of quantum chaos, and generally it was surmised 
14j that the stationary waves localized on one of these re- 
gions form an independent subset with specific dynamical 
properties. The frequency spectrum thus appears as the 
superposition of independent frequency subsets reflecting 
the phase space structure. This surmise has been found 
to be a good approximation for many systems, although 
some correlations may remain between the frequency sub- 
sets due to modes localized at the border between zones 
or due to the presence of partial barriers in phase space 

It is therefore important to know if this spectrum or- 
ganization is valid in the high frequency limit where it 
is supposed to hold, and even more important to assess 
if it is still relevant to the observable acoustic modes. 



We have thus numerically computed exact axisymmct- 
ric modes of Eq. (Q}, using the method described in [||, 
in the frequency range [u;i,12a;i], u>\ being the lowest 
acoustic frequency. To establish a link with the asymp- 
totic ray dynamics, we use a phase-space representation 
of the modes known as Husimi distribution [16[ . In order 
to compare a three-dimensional mode with the acoustic 
rays on a two-dimensional plane, the mode amplitude 
is first scaled by the square root of the distance to the 
rotation axis [7(. The Husimi distribution is then con- 
structed from a cut taken along the PSS: H(so,k ) — 
| / $'(s) exp(-(s - s )7(2A 2 )) exp(ifc s)ds| 2 . Here 
is the scaled version of the integral is taken along the 
curve r = r p and k is the moment in the direction tangent 
to this curve; A is the width of the Gaussian wavepacket 
on which is projected. In Fig. 2, (s 0l k Q ) is replaced 
by (9, kg) for comparison with data from Fig. 1. Ex- 
cept for some of the largest lengthscale modes close to 
the frequency wi, we find that the Husimi distribution 
enables to unambiguously associate the modes with the 
main structures of phase space. As illustrated in Fig. 
2, we distinguished the island modes trapped into the 
main stable islands, the chaotic modes localized in the 
central chaotic region and the whispering gallery modes 
associated with the whispering gallery region. We note 
that, due to the relatively low frequency considered, the 
chaotic modes do not spread over all parts of the chaotic 
region. 

Having defined subsets of modes, we can now analyze 
the properties of the corresponding frequency subsets. 
As shown in Fig. 3, the frequency spacings of the island 
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FIG. 3: Top: island mode frequencies in the interval 
[9oji, 12wi] showing the regular spacings corresponding to for- 
mula (5)j height reflects the value of the quantum number 
£. Bottom: Integrated spacing distribution for the chaotic 
modes in the same interval (the full line). Data correspond to 
around 200 modes from two symmetry classes. The dashed 
line is the result for the Gaussian Orthogonal Ensemble of 
Random Matrix Theory while the dotted line corresponds to 
the Poisson distribution characteristic of integrable systems. 



FIG. 2: (Color online) Comparison of numerically computed 
acoustic modes from |T} with phase space ray dynamics on 
the PSS. Three typical modes have been selected: island mode 
(blue/dark grey), chaotic mode (red/grey), whispering gallery 
mode (green/light grey). Left: spatial distribution of $' on 
one quarter of a meridian plane. Black lines are the nodal 
lines, full (resp. dashed) colored/grey lines are the level curves 
for positive (resp. negative) values. Magenta line at r — r p 
is the chosen PSS. Right top: cut of the island mode along 
the PSS, from the pole to the equator. Right bottom: level 
curves of the Husimi distributions of the three modes (for 
the same of value of A = 0.12r s (7r/2)), showing that they 
are mainly concentrated inside the main stable island (for the 
island mode), in the chaotic region (for the chaotic mode), or 
in the whispering gallery region (for the whispering gallery 
mode). 



modes display striking regularities which lead to the sim- 
ple empirical formula: 



u n i = nS n + £S( + a 



(2) 



where n and I are natural integers, S n and 5t are 
uniform frequency spacings. The a constant being fixed 
by a given island mode frequency, the formula proves 
sufficiently accurate to identify the other island modes 
among the whole set of computed frequencies. The phase 
space representation of these modes reveal that these 
regular patterns can be attributed to the existence of the 
stable island region in phase space. Although a zoom 
on this region would show a complex structure involving 
chaotic trajectories and chains of small islands nested 



between deformed surviving tori, these small scale details 
can be overlooked for the relatively large wavelengths 
considered here. To retrieve formula (J2]) and to find out 
how 5 n and Se relate to the properties of the star, we 
follow an approach inspired by the quantization of laser 
modes in cavities 17] . Indeed, far from the boundary, 



our problem can be translated into the propagation of 
light in an inhomogeneous medium, l/c s playing the 
role of the medium index. Close to the stable orbit, we 
can apply the paraxial approximation. In this case, it 
is known that the wave beam solution is [181 ] £) oc 
He(y/2^/w(a)) exp(-£ 2 /u;(cr) 2 ) exp(—i<f>(a, £)), where 
a, £ are coordinates parallel and transverse to the central 
periodic orbit (the yellow curve in Fig. 1), Hg is the Her- 
mite polynomial of degree £. The spreading of the beam 
in the transverse direction is described by w(a) which 
verifies (l/c s )d/da[(l/c s )(dw/da)) + a(a)w = 4/w 3 , 
where a(a) = (l/cl)d 2 c s /d£, 2 . The wave phase is 
ct>(a, 0) = w /„" da'/c s -2(£ + l) % c s da' /w 2 +£ 2 /(2c s i?) 
where R — w/(dw/da) is the radius of curvature of the 
beam wavefront. Numerically computed island modes 
have a transverse variation confirming this approxima- 
tion. We then obtain a stationary solution by imposing 
that the wave interferes constructively with itself. This 
requires that the phase accumulated following the 
periodic orbit (£ = 0) from one side of the boundary to 

the other side is ui J & do jc s — 2(£+ 1) J & c s da/w 2 = mr. 
This leads to the formula ^ with S n — it / {j^ da / c s ) 
and Si = 2(jk c s da/w 2 )/(j^ da/c s ). The numerical 
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value of S n obtained from the island mode frequencies 
shown on Fig. 3 (equal to 0.5514 in units of y/GM/r^ (0) 
where M is the stellar mass and G the gravitational 
constant) is well approximated, within 2.2 percent, by 
the theoretical one (equal to 0.5635 in the same units). 
While 5 n probes the sound velocity along the path of the 
periodic orbit, Si is obtained by solving the second order 
equation verified by w together with the two boundary 
conditions given by the necessity to match R with the 
radius of curvature of the two bounding surfaces. Thus, 
Se probes the second order transverse derivative of the 
sound velocity along the same path as well as the radius 
of curvature of the bounding surfaces. We note that 
similar modes around stable periodic orbit have been 
constructed in other systems, usually with the more 
systematic procedure of fin ding the normal forms and 
using EBK quantization li|, 19j . 

Having shown that modes whose Husimi distribution is 
localized in the near-integrable region display integrable- 
like quantization conditions, we now turn to the modes 
localized in the chaotic region. The subset of chaotic 
mode frequencies shows typical signatures of wave chaos 
such as frequency repulsion. Indeed, in Fig. 3, the in- 
tegrated distribution of consecutive frequency spacings 
Si = Wi+i — u>i (normalized by the mean frequency spac- 
ing of those modes) is much closer to the Random Ma- 
trix Theory result typical of chaotic systems [2(| than 
to the Poisson distribution result characteristic of inte- 
grable systems. This frequency statistics together with 
the fact that the corresponding modes are all localized 
in the chaotic region of the ray dynamics give a strong 
evidence that wave chaos occurs in rapidly rotating stars. 
The difficulty to solve Eq. ([I]) even with state of the art 
computational techniques prevents us to reach a larger 
frequency sample and to make detailed comparison with 
Random Matrix Theory as in e.g. (20j . 

The whispering gallery modes and the modes trapped 
in smaller island chains being associated with near- 
integrable regions of phase space, their frequencies are 
therefore expected to display regular patterns. The detail 
study of these regularities shall be considered elsewhere 
as it requires more mode calculations with a higher nu- 
merical resolution. It is also important to point out that 
these modes will be the most difficult to detect. Indeed, 
due to their small latitudinal wavelength (see Fig. 2), the 
positive and negative light intensity fluctuations strongly 
cancel out when integrated over the visible disk. 

Our results demonstrate that ray dynamics and quan- 
tum chaos concepts provide a qualitative and quantita- 
tive insight into the frequency spectrum of rapidly ro- 
tating stars. In particular we are able to separate the 
spectrum of a reasonably realistic star model into inte- 
grable and chaotic subsets. Being much less demanding 
than the direct eigenmodes computation as well as easily 
adaptable to non-axisymmetric modes (for which regular 
frequency patterns have also been found numerically 21|) 



and to more realistic stellar models, ray dynamics will be 
essential to further specify the asymptotic properties of 
the oscillation spectrum. 

The present analysis opens new perspectives in seis- 
mology of rapidly rotating stars. Observed spectra differ 
from theoretical ones as poorly understood mechanisms 
governing the intrinsic mode amplitude determine the 
frequencies that are actually detected. In this context, 
a priori information on the structure of the spectrum is 
crucial in order to identify the observed frequencies with 
specific stellar oscillation modes. Our results strongly 
suggest to first look for regular patterns to identify the 
island modes and to determine seismic observables as 8 n 
and 5g containing information on the star's interior. The 
remaining chaotic modes are also of special interest for 
seismology purposes: they are highly sensitive to small 
changes of the stellar model and, contrary to non- 
radial acoustic modes of slowly rotating stars which avoid 
the star's center, they probe this region which is crucial 
for stellar evolution theory. If enough chaotic modes are 
seen, their mean frequency spacing which is known to de- 
pend on the volume of the chaotic region would constrain 
the stellar rotation. 
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